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We report high-sensitivity AC susceptibility measurements of the penetration depth in the Meiss- 
ner state of the layered organic superconductor K-(BEDT-TTF)2Cu[N(CN)2]Br. We have studied 
nominally pure single crystals from the two different syntheses and employed controlled cooling 
procedures in order to minimize intrinsic remnant disorder at low temperatures associated with the 
glass transition, caused by ordering of the ethylene moieties in BEDT-TTF molecule at T G — 75 K. 
We find that the optimal cooling procedures (slow cooling of -0.2 K/h or annealing for 3 days in the 
region of T G ) needed to establish the ground state, depend critically on the sample origin indicat- 
ing different relaxation times of terminal ethylene groups. We show that, in the ground state, the 
behavior observed for nominally pure single crystals from both syntheses is consistent with uncon- 
ventional d-wave order parameter. The in-plane penetration depth X\ n (T) is strongly linear, whereas 
the out-of-plane component \ ou t(T) varies as T 2 . In contrast, the behavior of single crystals with 
long relaxation times observed after slow (-0.2 K/h) cooling is as expected for a d-wave supercon- 
ductor with impurities (i.e. A; n (T) oc A ou t(T) oc T 2 ) or might be also reasonably well described 
by the s-wave model. Our results might reconcile the contradictory findings previously reported by 
different authors. 



PACS numbers: 74.70.Kn, 74.25.Ha, 74.25.-q, 74.62. Bf 
I. INTRODUCTION 



Since the discovery of superconductivity in k-BEDT- 
TTF based materials a decade agcocl, the question of 
pairing symmetry has remained as one of the most in- 
triguing issues. From the very beginning, these mate- 
rials have attracted a lot of interest, not only because 
they achieved the highest superconducting (SC) tran- 
sition temperatures (Tc) among organic materials, but 
also because of their similarities to the high-temperature 
cuprate superconductors. First, the k-(BEDT-TTF) 2 A" 
(abbreviated as k-(ET) 2 A) are strongly anisotropic, 
quasi- two-dimensional materials, with a very weak in- 
terplane coupling. This feature is due to the crys- 
talographic structure, in which orthogonally aligned 
BEDT-TTF dimers form 2D conducting layers sand- 
wiched between the polymerized anion (X) layers. Sec- 
ond, antiferromagnetic (AF) and SC phases occur next 
to one another, which suggests that electron correla- 
tions play significant role in the establishment of the 
ground state. Indeed, the ground state of k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]C1 material is an insulating AF phase 
with mildly canted spinsEHa, while the ground state 
of K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br (abbreviated as k- 
(ET) 2 Br) and K-(BEDT-TTF) 2 Cu(NCS) 2 (abbreviated 
as k-(ET) 2 NCS) is a SC phase. In the former, the applied 



pressure supresses AF state and stabilizes SC stateS, 
whereas by deuterization of K-(ET) 2 Br, the ground state 
is gradually pushed from SC toward AF stateO. The 
phase diagram is, therefore, quite similar to that of the 
cuprates if doping is replaced by pressure or deuteriza- 
tion. Third, the normal state has some properties that 
are distinct from conventional metals, supporting the im- 
portance of the electron correlations. In particular, the 
Knight shift decreases significantly below about 50 K sug- 
gesting a suppression of density of states, that is, the 
appearance of a pseudogap near the Fermi energyQ. A 
broad dip in the electronic density of states around Fermi 
energy was also observed by,-STM measurements below 
about the same temperature^. Further, there is a peak 
in \jT{T at about 50 K, where T\ is the nuclear spin- 
lattice relaxation time, which suggests the presence of 
short range AF correlations!! This peak disappears un- 
der pressure concomitantly with superconductivity. 

The presence of significant electron correlations 
strongly favors the possibility of an unconventional SC. 
Results in favor of d-wave order parameter have been ob- 
tained by different experimental techniques. 13 C NMR 
measurements showed that the spin-lattice relaxation 
rate follows T 3 dependence at very low temperatures. 
This result, together with the Knight shift experiment, 
provide*, evidence for spin singlet pairing with-podes in 
the gapQ. The low temperature specific heatt°l, as well 
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as the thermal conductivity^, also showed a power-law 
behavior c s (T) oc T 2 and k(T) oc T, respectively. Fur- 
ther, magnetic penetration death A(T), measured by mi- 
crowave cavity perturba,tionE3, muon-spin relasationEj, 
tunnel diode oscillatorE-l and AC susceptibilitytSEJ, also 
displayed the power law behavior, usually in the form of 
a T and/or a T 2 behavior for A(T) at low temperatures. 

Recent angle resolved measurements of SC sap struc- 
ture using STMO and thermal conductivitytl clearly 
showed the fourfold symmetry in the angular variation, 
characteristic of the d-wave superconducting gap. Both 
measurements have revealed that nodes are directed 
along the directions rotated by 45 degrees relative to 
the in-plane crystal axes, indicating the d x 2_ y 2-wave su- 
perconductivity. Such a nodal structure indicates that 
both Fermi surfaces (oval-shape quasi-two-dimensional 
hole .cylinder band and an open quasi-one-dimensional 
bandt2l) should participate in SC pairing in contrast with 
theoretical predictions of superconductivity induced by 
AF spin fluctuations^. 

In jsoateast, some other penetration depth 
studiescilc30, as fXfA as the most recent specific 
heat measurementsOt3, led to results favoring con- 
ventional s-wave order parameter. In particular, a 
strong-coupling s-wave order parameter was claimed to 
be observed in the latter experiments. 

The question arises what is the source of the conflicting 
results and how could this discrepancy be resolved. As far 
as experimental determination of A(T) in the mixed state 
is concerned, the complex vortex dynamics might present 
a serious problem already at-fields as low as 70-300 Oe, 
as pointed out by Lee et alEfi. Further, additional com- 
plications might be due to an order-disorder transition 
that bears glassy features, taking place at Tq w 75 K 
for K-(ET) 2 B^|Or at T G1 w 70 K and T G2 w 53 K for 
K-(ET) 2 NC&iJ€a The transition region is situated be- 
tween 65 K and 85 K, and between 45 K and 75 K for 
K-(ET) 2 Br and k-(ET) 2 NCS systems, respectively. The 
transition is ascribed to the gradual freezing down of the 
remaining motion of the ethylene groups of the BEDT- 
TTF molecules that are thermally activated at high tem- 
peratures between the two possible conformations. That 
is, the relative orientation of the outer C-C bonds can be 
either eclipsed or staggered. Upon lowering the temper- 
ature, the former and latter are adopted for K-(ET) 2 Br 
and k-(ET) 2 NCS, respectively. X-ray diffraction mea- 
surements showed that at 125 K (iC|H 2 ) 2 groups are or- 
dered in average in the whole bulkEa. However, the pas- 
sage through the region of glassy transition appears to 
play a crucial role regarding the level of residual intrin- 
sic disorder at low temperatures. This might be also 
due to the anomalous changes in thermal exnaiision be- 
havior inside the same temperature regioriHEil. Rapid 
cooling rates are reflected in smaller resistivity ratio be- 
tween Tq ancLXis and larger resistivity humps centered 
at about 60 KQO The understanding of transport prop- 
erties in the normal state is further complicated by the 
fact that the standard resistivity behavior that resembles 



a semiconducting state above 100 KE3, whereas it be.- 
comes metallic below, is not reported for all syntheseso. 
In the latter case, samples display a metallic behavior in 
the whole temperature range between room temperature 
(denoted as RT) and Tq. In addition, indications are 
given that some Cu(II) may replace regular .Cu (I) during 
synthesis, affecting the resistivity behavior£3. The cor- 
relation between different mean free path in samples of 
different syntheses and presence of Cu(II) ions was also 
suggested!!! Further, sample dependence and relaxation 
effects were also observed iiHhe magnetization measure- 
ment results for K-(ET) 2 BrtZl. Samples of one synthesis 
show two peaks in the magnetization vs. field (M — H) 
curve, in contrast to crystals of another synthesis that 
show only one. Finally, the observed anomalous cooling 
rate dependence of M — H curve was attributed to the 
change in the resistivity curves and remnant disorder in 
the sample. 



In an attempt to reconcile the existing contradictions, 
and determine the pairing symmetry of the genuine SC 
ground state, we have undertaken an investigation that 
covered a broad range of single crystals of various syn- 
theses and in which the influence of thermal cycling and 
sample history was checked in carefully designed experi- 
ments. We have performed a high-resolution AC suscep- 
tibility measurements in the Meissner state in the two 
field geometries, i.e. when AC field was parallel and per- 
pendicular to the crystal planes. A quantitative data 
analysis, we have elaborated with scrutiny to account for 
the demagnetization correction in the latter geometry, 
enabled us to get a full characterization of each sample 
under study clarifying in this way a previously suspected 
sample dependence. Indeed, our experiments unfold that 
cooling procedures in the region of the glass transition, 
necessary to establish the bulk SC ground state, critically 
depend on the sample origin. In particular, our results 
demonstrate that the low temperature state is critically 
determined by time scale of experiment in the region of 
the glass transition, revealing in that way different relax- 
ation times of terminal ethylene groups in samples from 
two distinct syntheses studied in the most detail. We 
show that in the ground state, the behavior observed for 
nominally pure single crystals from both syntheses is con- 
sistent with unconventional <i-wave order parameter, that 
is the low-temperature in-plane and out-of-plane super- 
fluid density is proportional to T and T 2 , respectively. 
However, the behavior of the superfluid density observed 
after standard slow cooling for single crystals with long 
relaxation times deviates from that found in the ground 
state. It is as expected for a (i-wave superconductor with 
impurities or might be reasonably well described by the 
s-wave model. These results may give a solution to a 
decade long mystery in regard to the symmetry pairing 
in k-BEDT-TTF superconductors. A preliminary report 
was given in Ref. 
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FIG. 1: The K-(ET)2Br (above) and its reference sample (be- 
low). 



ensure that its dimensions and the shape were as close 
as possible to the original (see Fig. |l|). We start with 
the conventional formula, which relates the induced volt- 
age on the detection coil U, the measured susceptibility 
Xm; the demagnetization factor D, the volume V and the 
internal susceptibility x' , that is x' m corrected for the de- 
magnetization effect 1/(1 + Dx'): 



u^v x ' m = v 



X 



1 + D X ' 



(1) 



Here we have taken into account that measured samples 
are much smaller than the detection coil. If we use a con- 
ventional calibration assumption that the susceptibility 
for Niobium reference sample at low temperatures can 
be taken to be -1, we can finally determine the internal 
value of the susceptibility according to formula 



II. EXPERIMENTAL TECHNIQUE 

Measurements of the complex AC susceptibility (x — 
X' + were performed using a commercial AC suscep- 
tometer (CryoBIND/Sistemprojekt, Zagreb). The sen- 
sitivity of the system, expressed in equivalent magnetic 
moment, was Am = 2 • 10 -9 emu-in the broad tempera- 
ture range between 1.5 K and TcE?I Measurements were 
performed with H AG = 14 mOe at / = 231 Hz. Gold 
0.07% iron - copper thermocouple was used as a ther- 
mometer. The sample was placed in the upper one of 
the two identical secondary coils immersed in the liq- 
uid helium bath, positioned in such a way that applied 
AC field was either perpendicular or parallel to the con- 
ducting planes of the single crystals under study. In 
the remainder of the paper, the case of magnetic field 
aligned with the b crystallographic axis, i.e. perpendic- 
ular to the high conducting ac crystallographic planes, 
is denoted by Hag -L plane, while the case of magnetic 
field direction laying inside the ac crystallographic plane 
is denoted by H^c II plane. In order to probe the sam- 
ple in the Meissner state care was taken to reduce the 
amplitude of the AC field (Hag) until the component 
x'(T) was independent of Hac (H A g < 42 mOe) and 
the x"{T) component was negligible. No frequency de- 
pendence (13 Hz < / < 2 kHz) was observed for Hag < 
1 Oe. In addition no influence of the Earth's field was ob- 
served: runs performed with compensation of the Earth's 
field gave the same results. This is in accordance with 
the fact that the reported values for lower critical mag- 
netic field Hci(T), corrected for demagnetization, are 
far above the Earth's field He for all temperatures below 
8 K. In this temperature range Hci(T) > 10 Oecj, while 
the value of the Earth,'s field determined in our labora- 
tory is H E « 0.36 OeM 

The calibration of the system was performed with a 
piece of specially and carefully designed Niobium foil, 
which we will refer to, in the remainder of the paper, as 
the reference sample. For each K-(ET) 2 Br single crystal, 
a special reference sample was created independently to 



X = -77 



Vr 
V 



1-D t 



1+D X > 



u_ 



(2) 



where index r designates the reference sample. When 
the AC field was parallel to the crystal planes, the de- 
magnetization effect could be neglected, meaning that 
the middle factor in Eq. (||) can be simply replaced by 
unity. In order to correct the measured data for the de- 
magnetization factor for the AC field perpendicular to 
the crystal planes, we have developed a method, which 
strongly reduces the experimental error in x' to less than 
1% and allows us to get the reliable absolute value of the 
penetrated volume 1 + x' f° r this geometry, as well. The 
latter quantity directly determines the penetration depth 



and the s uperfluid density, as we show in Sections III B 



and |nC[ In the following, we argue and show that our 
claim is justified. 

Our starting point is that the reference sample repre- 
sents a perfect copy of the sample under study (Fig. [I]) , 
implying very close values of respective demagnetization 
factors D and D T . Since the K-(ET) 2 Br sample is almost 
completely diamagnetic in this geometry, we can use the 
approximation %' « — 1. Therefore, demagnetization ef- 
fects for both samples must be very close in value and 
the middle factor in Eq. (||) in the first order of approx- 
imation can be replaced by unity. Further, in order to 
improve precision and eliminate a small remaining dif- 
ference in demagnetization effect, the middle factor in 
the second order of approximation is calculated by the 
following procedure. Both K-(ET) 2 Br and the reference 
sample were taken to be fully superconducting disks (that 
is thin cylinders with aspect ratio of the length and the 
diameter of about 0.4) with the same face area and the 
same thickness as their originals. We consider the disk 
approximation to be more suitable for description of the 
real sample than the elipsoid one, used by authors of 
Ref. We base this assertion on the fact that the for- 
mer approximation better describes rather sharp sample 
edges, which might give a substantial contribution to the 
demagnetization factor. The middle factor in Eq. (||) is 
then given by the calculated ratio of diamagnetic effects 
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for these two bodies, 1 _p / prg~ ■ The systematic 
error due to the approximation of the specific shape of 
samples to the shape of the disk obviously cancels out by 
division. Numerical data for the demagnetization factor 
for the disk were taken from the literatureca. 



In order to calculate the area and the thickness, di- 
mensions of both K-(ET)2Br and reference sample were 
carefully measured with a high precision of 1%. The pre- 
cision was verified by the following procedure. All sig- 
nificant Niobium reference sample dimensions were mea- 
sured and the volume calculated. In addition, the same 
reference sample was weighted and the volume calculated 
using Niobium density, pNb = 8.57 g/cm 3 . The difference 
between two obtained values was always about 1%. Fi- 
nally, the middle factor in Eq. (^) obtained in this way 
differs from unity only for a few percent. This shows 
that the calculation procedure and the starting assump- 
tion are valid. 

In the end, we estimate that our calibration procedure 
for Hag -L plane is accurate in 1 + x' to 33% for x' close 
to -1, and to about 50% for -0.96 > x' > -0-5. As for 
the Hag II plane, accuracy in 1 + x' ls estimated to be 
about 15% for all measured low-temperature x! values. 

Nominally pure single crystals of K-(ET) 2 Br, origi- 
nating from two different syntheses had different re- 
sistivity fafjkis i?i?(75K/T m j n ) for similar cooling rates 
employedMEl. RR(75K/T min ) = p(75K)/ Pmin , where 
Pmin is the resistivity measured at temperature just above 
the SC transition. Single crystals with i?i?(75K/T m j n ) 
rs 200 and « 50 obtained by standard slow cooling rate 
have been labelled as SI and S2, respectively. 

Three different cooling procedures were used to cool 
samples from RT to 4.2 K. Special care was taken in 
the temperature rangeJaQ_K < T < 100 K, where relax- 
ation processes appearEZrE3. For Quenched (denoted as 
Q) state the sample was cooled down to liquid helium 
temperature in about one minute, which represents the 
average speed of about -300 K/min. over the whole tem- 
perature range. For Relaxed (denoted as R) state the 
sample was first cooled to 100 K in about 10 minutes. 
Between 60 K and 100 K the cooling rate was carefully 
monitored to amount to about -0.2 K/min. Below 60 K, 
the sample was finally cooled to 4.2 K in a few minutes. 
For the Annealed state (denoted as A) the sample was 
cooled down to liquid nitrogen temperatures in about one 
hour. Then it was kept between liquid nitrogen tempera- 
ture and 100 K for three days. Finally, it was cooled down 
to liquid helium temperature in a few minutes. Samples 
used in this study were rhombic platelets with face ar- 
eas between 0.51 and 2.15 mm 2 and thickness between 
0.29 and 0.70 mm. 



III. RESULTS AND ANALYSIS 



A. Complex susceptibility 



Our first important result concerns the influence of 
cooling rate on the components of the complex suscepti- 
bility in the SC state as a function of synthesis procedure. 
We present the behavior obtained in two principal field 
geometries. That is Hag -L plane and H^g || plane, as 
defined in Section |l[ On the basis of our previous mea- 
surements with the magnetic field aligned with the a and 
c crystallographic axesEj, we know that the specific ori- 
entation within the ac crystallographic plane does not 
influence the obtained result for the penetration depth 
since the out-of-plane component is a dominant contri- 
bution. Therefore, we always made sure that the field 
was aligned with the largest dimension of the platelet 
face in order to minimize the demagnetization factor. 

Susceptibility data obtained for sample SI for two dif- 
ferent cooling rates are shown in Fig. ||. We identify the 
R state, as the ground state. In the ground state, su- 
perconductivity sets-in at 12.0 K and 11.5 K, for Hag -L 
plane and Hag II plane, respectively. We have already 
reported difference in Tc, as well as anisotropy in the 
susceptibility for Hag perpendicular and parallel to the 
planealj. Here we report, to the best of our knowledge 
for the first time, the absolute susceptibility values in 
both field geometries. For Hac -L plane the sample re- 
sponse is almost completely diamagnetic, while for Hag II 
plane the susceptibility is somewhat smaller in magni- 
tude. Next, our new result points to a huge effect of 
cooling rate on the susceptibility value and Tc- When 
sample was cooled faster, the absolute value of suscep- 
tibility was smaller, and Tc lower. In other words, the 
diamagnetic region shrinks in x' vs - T plot. This huge 
effect is especially emphasized for Hag II plane geome- 
try, where the absolute x' value is almost an order of 
magnitude smaller for Q state than for R state. 

Susceptibility data obtained for sample S2 for three 
different cooling rates are presented in Fig. t| The first 
feature, reflecting a different sample quality, is that the 
ground state, characterized as before by almost complete 
diamagnetic response, is established in the A state, and 
not, as in the case of sample SI, in the R state. For a 
purpose of clarity, we will refex_.to the R state of sam- 
ple S2 as the intermediate stateLn. The anisotropy of the 
susceptibility in the ground state is somewhat larger for 
sample S2 than for sample SI. The second feature is in 
that the cooling rate effect on the susceptibility value is 
much smaller for sample S2 compared with the effects 
obtained for sample SI (see Fig. ^|). 

Finally, we point out that the x"(T) component was 
negligible, indicating clearly that the measured samples 
were in the Meissner state. 
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FIG. 2: Sample SI: Real and imaginary part of susceptibility 
for Relaxed (R) and Quenched (Q) state in (a) Hac -L plane 
and (b) Hac || plane geometry. The arrow in (a) illustrate 
the upper limit of the systematic error. 
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FIG. 3: Sample S2: Real and imaginary part of susceptibility 
for Annealed (A), Relaxed (R) and Quenched (Q) state in 
(a) Hac -L plane and (b) Hac || plane geometry. The arrows 
in (a) illustrate the upper limit of the systematic error. 
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Sample S2, H AC \\ plane (b) 








B. Penetration depth 



Further, we analyse susceptibility data for hereabove 
described single crystals of different quality in order to 
get the penetration depth temperature behavior and the 
values at zero temperature. In Hac II plane geometry we 
have already argued that the in-plane penetration depth 
Ai n can be neglected in analysis and that the out-of-plane 
penetration depth A out can be obtained from the suscep- 
tibility data using the formula for a thin superconducting 



plate in a parallel fieldt-3 (see Fig. £|(a)) 

1 + X'=^tanh(^) (3) 

'out V^out/ 

Now we address in more detail 7?ac -L plane geometry. 
The magnetic field is strictly perpendicular to the con- 
ducting planes, so is the responding magnetization that 
expels it out of the bulk. The resultant circulating su- 
percurrents will therefore flow within planes, which will 
only give contribution to Ai n (see Fig. 0(b)). In order to 
obtain an appropriate formula for the analysis, we start 
with the generalization of Eq. (||) . We take into account 
that the ratio between the doubled penetration depth 
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FIG. 4: Simplified depiction of the sample in the magnetic 
field (a) H\c || plane and (b) H\c -L plane. Shaded parts of 
the sample represent the volume penetrated by the magnetic 
field. Condition A out /A in 3> UirHout in (a) ensures that Ai n 
can be neglected in the analysis^. 



2A ou t and the sample width in the direction of field pen- 
etration / ut is the ratio between the volume penetrated 
by the magnetic field Vp and the whole sample volume 
V. Therefore, 



l + X '= V tanh-. 



(4) 



We note that, if the sample has the shape of the 
platelet, the magnetic field penetrates along the whole 
edge of the face. We can rewrite Eq. (0) as 



1 + x = — : — tanh ■ 



.4 



CA i; 



(5) 



where C is circumference, A area of the platelet face and 
Ai n the in-plane penetration depth (see Fig. ||(b)). 

Our second important result concerns the tempera- 
ture dependence of Ai n and A ou t in the ground state as 
a function of the synthesis procedure. In Figs. || and || 
we show the influence of cooling rate on the anisotropic 
penetration depth as a function of the sample quality. 
Note that in the ground state of both samples SI and S2 
the temperature dependence of Aj n and A ou t at tem- 
peratures below about 5 K is well described by the T 
and T 2 law, respectively. The upper bound of the fit 
range (5 K) is given by the general requirement that 
the genuine low-temperature behaviour of the penetra- 
tion depth is strictly obeyed only far from the critical 
region close to Tq. The full lines correspond to the cal- 
culated fit to the power law behavior in the temperature 
range 1.6 K < T < 5 K 



A,-, 



/••in I £ 



T 



Ain(O) 



fcout ( 7fT ) +Aout(0). 



(6) 
(7) 



We get k in = 5.2 /im, A; n (0) = 1.5±0.5/an, k out — 56 /im, 
Aout(O) = 53 ± 10/zm and k in = 2.8 fim, A in (0) = 1.1 ± 
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FIG. 5: Sample SI: (a) In-plane and (b) out-of-plane pene- 
tration depth for Relaxed (R) and Quenched (Q) state. Full 
lines represent the fit to the power law behavior, while the 
arrow in (a) illustrate the upper limit of the systematic error. 



0.4/xm, fc out = 69 /im, A out (0) = 85±10^mfor SI and S2, 
respectively. The penetration depth values at K, A(0), 
observed in the ground states of both samples, are in 
a very good accordance with^alj^|4tp,jia-j>«netration 
depths given in the Uterate^yMakUdMa. 

A special attention should be given to the relative 
change of the penetration depth at low temperatures, 
r}(T) = (A(T) - A(0))/A(0). We denote the deviation of 
Aj n (T) and A out (T') from their values at K, in units 
of A(0), as i]- m (T) and rj out (T), respectively. We find 
that J&a(5K) = 1.4, r? out (5K) = 0.19 and ?y in (5K) = 1.1, 
r? 0Ut (5K) = 0.15 in the ground state of samples SI and S2, 
respectively. First we note that r\ values are in perfect ac- 
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FIG. 6: Sample S2: (a) In-plane and (b) out-of-plane pene- 
tration depth for Annealed (A), Relaxed (R) and Quenched 
(Q) state. Full lines represent the fit to the power law behav- 
ior, while the arrows in (a) illustrate the upper limit of the 
systematic error. 



cordance for both samples confirming that they are in the 
same - ground - state. Second, a difference between if m 
and rj out values for almost an order of magnitude proves 
a strong anisotropy in physical properties between the 
two orientations. This also confirms our choice for the 
exponent in the power law in Eqs. (||) and ((?]). 

Finally, our third important result concerns the inter- 
mediate state (R state in sample S2), for which we get 
77 in (5K) = 0.13, ??out(5K) = 0.10. Unlike for the ground 
state, the temperature dependences of A; n and A out are 
so similar that they cannot be any longer described by 
different power laws as was the case in the ground state. 



The obvious solution is to try to fit both penetration 
depths to the T 2 power law behavior 



( T V 

Aout = fcout I 7jT J + Aout (0) i 



(8) 
(9) 



which suggests a G?-wave superconductor with impuri- 
ties. We get k- m = 9.2 fim, Ai n (0) = 12 ± 6^m, and 
feout = 58 jum, Aout(0) = 110 ± 20^tm. We point out 
that the fit to the s-wave model describe our data al- 



most equally well (see Section III C ). On the other hand, 
the s-wave model fails completely for the penetration 
depth temperature dependences in the ground state of 
both samples SI and S2. 

Now we comment penetration depth results for Q state. 
We find that r; in (5K) = 0.11, % u t(5K) = 0.12 and 
77in(5K) = 0.16, 77out(5K) = 0.10 for samples SI and S2, 
respectively. Here we apply the same arguments as in 
the case of the intermediate state and fit both Ai n and 
A ou t to the T 2 power law behavior (Eqs. (||) and We 
get A in (0) = 100 ± 50^m, A out (0) = 830 ± lOO^rn, and 
A in (0) = 24± 12/im, A out (0) = 170± 20^m for sample SI 
and S2, respectively. The fact that these fits describe 
the penetration depth data well again suggests a rf-wave 
superconductor with impurities. Further, it should be 
noted that (i) A(0) values are larger for Q state than 
for the intermediate state, suggesting a larger disorder in 
the former state, (ii) A(0) values for sample SI are signif- 
icantly larger than ones for S2, suggesting significantly 
larger disorder in the former sample in Q state. Finally, 
the result that A ou t(0) for sample SI is close to the crys- 
tal size indicates that the bulk superconductivity is not 
established, allowing us to define the boundary between 
bulk and nonbulk SC at x' = —0.7 for Hag -L plane. 



State 


Ain(0) 


Ao Ut (0) 




&out 


J?in(5K) 


r?out(5K) 




(jun) 


(jan) 


(fim) 


(fim) 






Sample SI 


R 


1.5 ±0.5 


53 ± 10 


5.2 


56 


1.4 


0.19 


Q 


100 ± 50 


830 ± 100 


60 


480 


0.11 


0.12 


Sample S2 


A 


1.1 ±0.4 


85 ± 10 


2.8 


69 


1.1 


0.15 


R 


12 ±6 


110 ±20 


9.2 


58 


0.13 


0.10 


Q 


24 ± 12 


170 ± 20 


22 


88 


0.16 


0.10 



TABLE I: Penetration depth properties, as defined in Text, 
for (i) Sample SI in Relaxed (ground) and Quenched state, 
(ii) Sample S2 in Annealed (ground), Relaxed (intermediate) 
and Quenched state. 



C. Superfluid density 

In the following, we address the temperature depen- 
dance of the superfluid density in order to get informa- 
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FIG. 7: Sample SI: (a) In-plane and (b) out-of-plane super- 
fluid density for the ground (R) state. Solid line is a fit to the 
polinomial expression. A large noise in p Sj ; n is due to small 
values of A; n (0) (see Text). 



FIG. 8: Sample S2: (a) In-plane and (b) out-of-plane super- 
fluid density for the ground (A) state. Solid line is a fit to the 
polinomial expression. A large noise in p Si ; n is due to small 
values of A; n (0) (see Text). 



tion on the symmetry of the superconducting state. We 
construct the in-plane superfluid density p Sj j n and out-of- 
plane superfluid density p s , ut as 



Ps,out 



Ain(O) 
Ain(T) 

Aout(O) 
Aout(T) 



(10) 

(11) 



p Sj i n and p Si out for the ground state of sample SI (es- 
tablished in the R state) and for the ground state of 



sample S2 (established in the A state) as a function of 
reduced temperature t = T /Tq are displayed in Figs. 
and ||, respectively. There is a strong resemblance in the 
behavior found for both samples. Note that the leading 
terms, which describe the low temperature behavior, are 
T and T 2 term for p Sj j n and p s ,out ; respectivelyO. This is 
to be expected, because T and T 2 terms describe the low 
temperature behavior of the A; n and A out in the ground 
state. If we expand the leading term to the full poly- 
nomial, in order to fit the superfluid density data in the 
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whole temperature region below Tq , we finally obtain for 
sample SI 

p s , in = 1 - 1.95t + 1.45i 3 - 0.09i 4 - 0.41i 5 (12) 
p s out = 1 - l.88t 2 - 0.73i 3 + 4.47t 4 - 2.86i 5 , (13) 

and for sample S2 

Ps.in = 1 - 1.68* + 0.78i 3 + 0.16i 4 - 0.26t 5 (14) 
Ps.out = 1 - l-45i 2 + 2.98i 3 - 5.38i 4 + 2.84i 5 . (15) 

Taking into account a relatively large experimental er- 
ror in the penetration depth values (see Figs. || and ||), 
the leading coefficient values might be considered to be 
almost the same. In addition to the systematic error, 
there is also a noise, which is a mere consequence of the 
fact that p Sj i n is calculated according to the expression 
Ps,in(T) = (Ai n (0)/Aj n (T)) 2 , so that the absolute noise in 
/o Sj i n is proportional to the relative noise in the penetra- 
tion depth data. That implies a larger noise for smaller 
values of Ai n (0), which becomes substantial for Aj n (0) of 
the order of lfim. Finally, we point out that the shape of 
curves for the in-plane and the out-of-plane penetration 
depth in the ground state are qualitatively different from 
the s-wave dependence. 

The behavior for the in-plane p Si j n and out-of-plane su- 
perfluid density p Si0u t for the intermediate state of sam- 
ple S2 (established in R state) is displayed in Fig. ||. 
Here, p Sj i n data are insensitive to the systematic error at 
low temperatures, since the correction for Ai n drops out 
from the expression for p Sj i n . Note that the leading term, 
describing the low temperature behavior, is T 2 term for 
both the in-plane and out-of-plane superfluid density: 

p s in = 1 - 1.69i 2 + 2.77i 3 - 5.05i 4 + 2.97i 5 (16) 
p s>out = 1 - 1.02t 2 + 1.34i 3 - 2.50i 4 + 1.18t 5 (17) 

In order to demonstrate the fact that the superfluid den- 
sity behavior in the intermediate state is also rather close 
to the dependence expected for the s-wave order param- 
eter behavior, the s-wave model dependence is added as 
a dashed line for both orientations. 



IV. DISCUSSION 

We start discussion by pointing out that the well- 
defined, ground state properties - complex susceptibil- 
ity, penetration depth and superfluid density - were es- 
sentially reproducible for all studied single crystals from 
both syntheses SI and S2. Cooling rate dependent ef- 
fects were also reproducible, but the observed behavior 
was the same only for single crystals from the same syn- 
thesis, while it differed significantly from the observed 
behavior for single crystals from the other synthesis. 

First, we would like to comment on the appotropy in 
T c . As in the previously published papertiil, we have 
established that ATc = 0.5 K for Hac -L an d -Hac II 




FIG. 9: Sample S2: (a) In-plane and (b) out-of-plane super- 
fluid density for the intermediate (R) state. Solid line is a 
fit to the polinomial expression, and dashed line presents the 
s-wave model (see Text). 



plane geometry cannot be ascribed to the experimental 
error. We suggest that this anisotropy might be the con- 
sequence of the fact that the diamagnetic shielding is no 
longer effective for Hac II plane geometry in the range of 
0.5 K below Tc, which is due to small sample dimensions 
and huge out-of-plane ( Josephson-like) penetration depth 
near Tc. However, we point out that Tc anisotropy has 
negligible effects, if any, to our data analysis, which is pri- 
marily done in the low-temperature region, and therefore 
does not influence the resulting conclusions. 
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Second, we address the observed differences in cool- 
ing rate effects between samples SI and S2. As pointed 
out earlier, both samples show the same behavior in the 
ground state: almost full diamagnetism for -Hac -L plane, 
the same temperature dependence and zero-temperature 
value of the in-plane and the out-of-plane penetration 
depth, as well as the same superfluid density tempera- 
ture dependence. However, we note that sample S2 re- 
quired a completely different cooling procedure with sig- 
nificantly longer time spent in the temperature region 
around 80 K, compared to sample SI, in order to ac- 
complish the ground state. Moreover, sample S2 is much 
less sensitive to the cooling rate. That is, the differ- 
ence regarding low-temperature susceptibility and zero- 
temperature penetration depth values between A and 
Q states for sample S2 is much smaller than the difference 
between R and Q states for sample SI. Both features 
indicate that (i) the low-temperature state is critically 
determined by the time scale of experiment in the region 
of glass transition and (ii) relaxation times of ethylene 
groups in the single crystals originating from synthesis S2 
are much longer than the ones in the single crystals orig- 
inating from synthesis SI. When the applied time scale 
is much longer than the relaxation time of the ethylene 
moieties, the low-temperature state is the ground state. 
In contrast, if the relaxation time exceeds the time scale 
of experiment, remnant disorder at low temperatures will 
be substantial and Q state will be established. Different 
relaxation times of ethylene groups might also explain 
why the resistivity ratio i?i?(75K/T m i n ) is much larger 
for samples from synthesis SI when standard slow cooling 
rate is applied. At this stage, we can only speculate about 
the possible origin of different relaxation times. The ex- 
perimental observations that the crystals from synthe- 
sis S2 show weak metallic behavior, instead of a semi- 
conducting behavior between RT and 100 K observed 
for samples of synthesis SI, might be of the same origin. 
Since the RT resistivity values do not differ substantially, 
we propose that the subtle local variations of impurity 
level in nominally pure samples from different syntheses 
might be responsible for the observed differences. 

Next we comment on the behavior of the in-plane su- 
perfluid density in the ground state. The temperature de- 
pendence of the in-plane superfluid density for the d-wave 
superconducting order parameter A(k) = Acos(20), 
where 4> ls the angle between the quasiparticle momen- 
tum k aad,rkhe a axis, within the weak coupling theory is 
given byt 5 " 



1 - 0.64784- 0.276^ 



(18) 



The coefficient a of the leading term 4 in p Sj j n = 1— at+. . . 
depends strongly on the ratio of the superconducting 
transition temperature and the zero temperature super- 
conducting order parameter. A comparison of values for 
a in Eq. ( |l8|) to those in Eqs. ( |i~2|) and (14) suggests that 
the superconducting order parameter at T = K is much 
smaller than that predicted by the weak-coupling limit. 
As a result, this would also imply that the nodal region, 



a 

in 

a. 




T/T n 



FIG. 10: Superfluid density in d + s-wave model for a few 
values of the s-wave component parameter r. 



which volume is inversely proportional to the angular 
slope of the gap near the nodeea [i = 1/ A-dA(4>)/d~4>\node, 
occupies a much larger fraction of the phase space at low 
temperatures. We point out that we have reached the 
same conclusion already on the basis of previously ob- 
tained resultsO. In the latter case, we did not evalu- 
ate Ai n (0) from our measurements, but rather we used 
the values reported in the literature. Following the same 
procedure as in Ref. [l5|, we can again use uncalibrated 
data from present measurements and calculate the devi- 
ation of A; n from the minimum value at-the lowest at- 
tainable temperature A; n (T) — Ai n (T m i n )lij. Then, if we 
take the .penetration depth value at K given in the 
literatureE2loc§c3, Ai n (0) ~ 1/xm, we get the same re- 
sult, that is the leading coefficient is much larger than ex- 
pected in the weak-coupling model. In addition, we point 
out that the same behavior ofpOin for Aj n (0) < 1.3/Ltm was 
reported by Carrington et al&3. They pointed out that 
only at A; n (0) > 1.8/xm the slope becomes similar to the 
one reported for high-Tc cuprate superconductors and 
expected in the weak coupling model. From our present 
data, we get A; n (0) « 3/um for the crossover in-plane pen- 
etration depth value. 

One possibility to interpret our results is to con- 
sider the mixture of d-wave and s-wave order parameter, 
which corresponds to the superconducting order parame- 
ter A(fc) =rA(cos(20)+r), with r representing the s-wave 
componentcHl. In this case the leading linear coefficient a 
increases with the increase of r according to the expres- 
sion: 



2 In 2 



2.14V1 - r 2 



exp 



2r 2 



l + 2r 2 



(19) 



The shapes of the superfluid density curves for several 
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values of parameter r are given in Fig. [lfl. For our re- 
sults, \r\ 0.7 gives a very good agreement, which is, on 
the other hand, theoretically very unlikely. Recently, an 
admixture of s-wave component with —0.067 for k- 
(BEDT-TTF)Cu(NCS) 2 was suggested^ from the anal- 
ysis of thej-angular dependent magnetothermal conduc- 
tivity dataEJ. These data suggest that the d x 2_ y 2-wave 
instead of d xy -wave pairing is responsible for supercon- 
ductivity in the re-(ET) 2 A materials. This implies that 
both Fermi surfaces .participate in pairing in contrast to 
previous assumption^ Since in the former case the spin 
fluctuation model for the k-(ET) 2 A materials no longer 
describes the superconductivity pairing, the origin should 
be found elsewhere. Coulomb interaction, which is re- 
sponsible for the d-wave superconductivity, gives rise to 
both spin and charge fluctuations, so the obvious solu- 
tion appears to be that charge fluctuations play the prin- 
cipal role in the k-(ET) 2 X superconductivity. The value 
r = —0.07 suggests that the nodal lines in A(k) pass 
through the gap between two Fermi surfaces. This is 
consistent with the d + s-wave model in which the super- 
conductivity is due to charge fluctuations between differ- 
ent groups of dimers. On the other hand, for r m —0.7, 
the nodal directions cross the 2D circular Fermi surface, 
and for r w 0.7, the nodal directions cross the ID Fermi 
surface. If d + s superconductivity model is generated by 
charge fluctuations, such a scenario is unlikely to work, 
since this implies the strong intra Coulomb repulsion in 
each energy band. 

Now, we address the behavior of the in-plane super- 
fluid density in the intermediate state (Fig. ||). Con- 
sidering the polynomial fit in Eq. (|l6|), we see that the 
in-plane superfluid density appears to fit very well the 
d-w»¥«[-m©del with impurity scattering in the unitary 
limitEl'Eal^l. However, there is a serious discrepancy be- 
tween our experimental data and the impurity model. 
First of all, we ascribe the difference between A state 
(ground state), R state (intermediate state) and Q state 
to the residual degree of the ethylene disorder. Annealed 
state we discussed is the ground state with the least dis- 
order in the ethylene groups. It is natural to assume that 
the ethylene disorder gives rise to the quasiparticle scat- 
tering, which changes the i-linear dependence of p Sj j n into 
the i-squared dependence. But this is in contradiction 
with the well establishedrfact that a small disorder de- 
presses Tq dramaticallyEM3. In our experiment, it seems 
that Tc is practically unaffected by the ethylene disor- 
der for the intermediate state, achieved by slow cooling 
of -0.2 K/min. in the region of the glass transition. Fur- 
ther, taking into account the fact that the A(0) increases 
for at least by a factor of 6, superconducting electron 
density at K, n s (0) cx A _2 (0), decreases at least by a 
factor of 36. A simple impurity model cannot describe 
the surprising combination of these two features. On the 
other hand, the results for other low-temperature states 
achieved by cooling rates qc < —1 K/min. for sample SI, 
and for Q state in sample S2 are more consistent with 
the theory. In these cases both Tc, x'(0) and therefore 



A(0) are concomitantly influenced by the remnant ethy- 
lene disorder. However, it is still difficult to correlate the 
observed behavior to the impurity model quantitatively. 
This discrepancy might indicate that the degree of disor- 
der at low temperatures, as defined by the cooling rate 
in the region of the glass transition, has also a profound 
influence on electronic correlations, responsible for SC 
pairing. At this point we would like to recall the result 
of Kund et alxB showing changes in the crystal structure 
parameters in the region of the glass transition. It might 
be that these changes are also susceptible to the cooling 
rate. In addition, authors of Ref. ^ have pointed out that 
the role of disorder in this class of superconductors might 
be different than in the other unconventional supercon- 
ductors due to the vicinity of the Mott insulator in the 
phase diagram. 

Coming back to the intermediate state, we suggest 
that the indication of the decrease in the superconduct- 
ing electron density may be related to the reduction of 
the superconductivity volume. It has been already re- 
ported that cooling rates combined by progressive deuter- 
ation influence the .taw-temperature electronic state in k- 
(ET) 2 Br samplesB'Eil. Deuterated K-(ET) 2 Br system is 
situated in the critical region between an insulating AF 
transition at 15 K and a SC transition at 11.5 K. Despite 
the slow cooling rate, the deuterated sample gives almost 
the same Tc as in the hydrogenated K-(ET) 2 Br system, 
the SC state is not fully established in the bulk. Note 
that this results strongly resembles what we observed for 
S2 samples in the intermediate (R) state. In addition, 
a gradual decrease of the susceptibility below Tc in the 
deuterated system strongly indicates the inhomogenous 
nature of the SC state. This is in contrast to what we ob- 
serve in sample S2, in which the susceptibility curves are 
rather sharp even in Q state. More rapid cooling rates 
induce a decrease of Tc and a substantial decrease of SC 
volume fraction. Authors of Refs. |(||54] have argued that 
since the electronic specific heat of rapidly cooled deuter- 
ated samples did not show any finite electronic contribu- 
tion at low temperatures, the missing part of supercon- 
ducting phase should be considered to be the magnetic 
insulating phase. The question arises if their conclusion 
might also be valid in the hydrogenated system. Tak- 
ing into account the structure of the phase diagram of 
this class of superconductors^, we think that this might 
be the case. However, specific heat data under carefully 
controlled cooling cycles are needed to resolve this issue. 

The most intriguing fact about Fig. [9] is that in the 
intermediate state the observed data could be well de- 
scribed by the s-wave model as well. This gives a possible 
explanation for the contradictory findings in favor of the 
s-wave and d-wave model in the same material. The be- 
havior is obviously strongly influenced by both thermal 
history as well as synthesis, which suggests that the same 
material was not measured in the same low-temperature 
state. We hope that our results could contribute to rec- 
oncile contradictory findings met in the past. 
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V. CONCLUSION 

The level of residual disorder and electronic properties 
at low temperatures are critically determined by the time 
scale of experiment in the region of the glassy transition 
and the sample synthesis. This fact imposes an addi- 
tional requirement to get a reliable description of the SC 
state and that is to perform a full characterization of 
the SC state in the sample under study in the same well 
defined and controlled cooling conditions. The origin of 
the observed differences can be attributed to the residual 
ethylene disorder, which might be theoretically consid- 
ered as the impurity effect. 

The in-plane superfluid density of the ground state 
with the lowest residual ethylene disorder exhibits clear 
T-linear dependence, which is consistent with the d-wave 
model and in contradiction with the s-wave model. The 
leading T-linear coefficient is much larger than the one 
expected for the weak-coupling d-wave model. The d+s- 
wave model can remove this numerical discrepancy. How- 
ever, rather large s-wave component needed to fit our 
data is not consistent with the recent thermal conductiv- 
ity results. 

On the other hand, the in-plane superfluid density 
of the intermediate state achieved by slow cooling of - 
0.2 K/min. in samples of one synthesis, as well as in the 
states achieved by rapid cooling rates (qc < — 1 K/min.) 
in all studied samples, clearly exhibits the T-squared de- 
pendence, consistent with d-wave superconductivity in 



presence of impurities. In the intermediate state, the 
same Tc as in the ground state is accompanied by a rel- 
atively large reduction in the superconducting electron 
density, which cannot be explained within the standard 
impurity model. In the rapidly cooled states, Tc does re- 
duce concomitantly with A(0), still the depression of T c 
is quantitatively too small and the reduction of the su- 
perconducting electron density is too large. More work, 
both theoretical and experimental, is needed to be done 
to resolve these issues. 

Finally, the in-plane superfluid density data in the in- 
termediate state can be relatively well fitted to the s-wave 
model data as well. This fact gives a possible explanation 
for the contradictory findings in favor of s-wave and d- 
wave model in the same material and will hopefully con- 
tribute to a long awaiting consensus regarding the pairing 
symmetry in the k-(BEDT-TTF) based superconductors. 
In this circumstance, specific heat measurements in the 
same well defined and controlled cooling conditions are 
highly desirable. 
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